Self-consistent self-energy analysis of photoemission data 
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Here we present the details of a self-consistent procedure of the photoemission data analysis 
within the self-energy approach introduced in Ref . (cond-mat/0405696). We derive the relations 
of the quasiparticle self-energy with the parameters provided by photoemission spectra, demonstrate 
self-consistency of the Kramers-Kronig procedure and its robustness in determination of the bare 
dispersion, examine the possible influence of a particle-hole asymmetry, discuss the necessity of a 
clear definition of the "kink", and demonstrate the applicability of the developed approach for a 
couple of samples. 
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INTRODUCTION 

Considerably improved characteristics of the electron 
analysers used in modern angle-resolved photoemission 
spectroscopy (ARPES) together with the accumulated 
experience in photoemission on superconducting cuprates 
as well as recent success in clarification of the under- 
lying electronic structure ^ H, II] have led to under- 
standing that a careful analysis of the lineshape of pho- 
toemission spectra 0: Hi is needed to clarify the origin of 
the main interaction which is responsible for the pairing 
mechanism. The nodal direction in the Brillouin zone 
(BZ) of superconducting cuprates, being responsible for 
their conducting properties, is of a special importance. 
While the pronounced doping dependence of the quasi- 
particle spectral weight in the antinodal region of the BZ 
unambiguously points to a magnetic origin of the strong 
electron-boson coupling seen by ARPES Hill [13, the 
nodal direction, where the famous "kink" on the renor- 
malized dispersion ll^l Q| has been shown to per- 
sists in the whole doping range 0, H^. is a subject of 
unresolved controversy between phonon [iSlllij and spin- 
fluctuation ^3 scenarios. 

Recently we have shown that the spectral weight 
of electronic excitations detected by ARPES along the 
nodal direction can be self-consistently described within 
the quasiparticle self-energy approach 0. With the sug- 
gested procedure the real and imaginary parts of the self- 
energy in a wide frequency range as well as the bare band 
dispersion can be uniquely determined. We have shown 
that the procedure is applicable to an underdoped Bi- 
2212 in the pseudo-gap state and determined a number 
of correlation parameters for that sample. 

Here we give a detailed description of the procedure, 
discussing its limitations and physical consequences. We 
derive the relations of the quasiparticle self-energy with 
the parameters of photoemission spectra, demonstrate 
self-consistency of the procedure and its robustness in 
determination of the bare band dispersion, examine the 



possible influence of a particle-hole asymmetry, discuss a 
necessity of clear definition of the "kink", and demon- 
strate the applicability of the developed approach for 
other samples: an overdoped Bi(Pb)-2212 and optimally 
doped Bi(La)-2201. 



SELF-ENERGY APPROACH 

We start with a definition of the spectral function 
within the self-energy approach: 
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Here e(k) is the bare band dispersion and S = E' + iE" is 
a quasiparticle self-energy, an analytical function the real 
and imaginary parts of which are related by the Kramers- 
Kronig (KK) transformation T^: 
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where PV denotes the Cauchy principal value. Within 
such a definition, T,"{ui) < 0, and S'(a;) > for w < 0. 

It is instructive to express some interaction parameters 
via both self-energy functions. For example, the coupling 
strength A, which can be defined as 
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can be expressed in terms of S" differentiating the KK 
relation 0: 



dj:'{uj) 

duj 



TT 



Y."{x)-^"{uj) 

{x — 



dx. 



(5) 



Here we used the fact that adding some constant to 

in (jSJ does not change the result. Then, for an even 
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FIG. 1: Real (blue) and imaginary (red) parts of the self-energy related by Kramers-Kronig (KK) transform: E' — KKE", 
for three models of E" tails. 
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Eq. ® gives a certain feeling how S"(a;) function in- 
fluences A. For example, for a simple case 
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where > is an energy cut-off and C = — S"(0) > 
is an offset, Eq. © gives 
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for C ojc- As we show below, the different but rea- 
sonable tails of (at \uj\ > ujc) almost do not effect 
the experimentally determined A/a ratio, i.e. influence of 
the high energy region on A can be described by only one 
parameter, Uc- 



NODAL SPECTRA ANALYSIS 

The self-energy approach can be app lied to the nodal 
spectra analysis (e.g., see Refs.|il[ll[ll as far as the 
seed dispersion in ^ along the nodal direction is iden- 
tical to be the bare electron dispersion, not affected by 
either superconducting gap or pseudo-gap opening. 

Modern ARPES measures the photocurrent intensity 
of outgoing electrons, /, which is the quasiparticle spec- 
tral weight (the quasiparticle spectral function A multi- 
plied by the Fermi-function /) slightly modified by ma- 
trix elements, M(k,hi'), energy and momentum resolu- 
tions, Roj{lo) and i?fc(k), and an extrinsic backgroimd 

/(w,k) (X MA{Lj,k)f{Lj) ® ^ Rk + B{lj). (9) 



It is peculiar for all cuprates, and especially for Bi- 
cuprates, that the matrix elements variation with k along 
the nodal direction is negligib le and that the extrinsic 
background is fc- independent '2C| , thus can be easily sub- 
tracted. The energy and momentum resolutions produce 
rather small effect which we discuss later. So, the pho- 
tocurrent intensity with subtracted background is simply 
proportional to the quasiparticle spectral weight, and the 
question is whether it can be completely described by 
Eq. CJ. 

The self-energy in Eq. |^ is believed to be momentum 
independent along the nodal direction, or at least much 
less dependent than on frequency: dE/cZfc <C VFdTj/duj, 
where vp = {de/dk)k=kp is the bare Fermi velocity. 
This assumption is supported by a perfect Lorentzian 
lineshape of the momentum distribution curves (MDC's, 
MDC(fc) = A{k)^^const) [3, although it has been 
pointed out [i^ that a linear fc-dependence of S' is a more 
general case which leads to the Lorentzian lineshape of 
MDC's. We will discuss a possibile fc-dependence of E' 
below. 

As long as MDC's are Lorentzians, two quantities at 
each frequency can be derived from the experimental 
data: the MDC width, W{uj) (defined as the half width 
at half maximum), and the MDC peak position, fc„i(a;), 
which is thus an inverse function to the renormalized 
dispersion. In the vicinity of the Fermi level one can 
consider the bare dispersion as linear, and the relations 
between the self-energy parts and these experimentally 
determined functions are especially simple 

I]'(a;) — Lo ~ vp[km(y>) — kp] 
I]"(w) = -vfW{lj). 
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Since vp is a. priori unknown, the system (|10|I - H11|) is 
incomplete. One can complete the system by Eq. ^ 
but encountering a problem of " tails" , the high energy 
behavior of that is also unknown. For example. 
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FIG. 2: Illustration of the fitting procedure: real parts of the self-energy T.'^i^p{u!) (filled blue squares) obtained by 1131 and 
E^jf (tj) (open violet circles) by H14^ and Q; the difference AT,' (to) — E^^(a;) — E^j^(cj) (small crosses) is fitted to R'{ui) 
(corresponding solid line), the contribution of overall resolution determined by 1171 - 118^ . In first three panels ujq = —0.9 but 
different n = 3, 4 and 6 in l|15|l are compensated by different uic = 0.34, 0.45 and 0.52 eV respectively. Last two panels, the 
"best" fitting results for slightly different ljo's. 



using another ultimate model for S"-tails, 



SELF-CONSISTENT KK-TRANSFORM 



S"(c.) = - 



C for \uj\ < ijjc, 
C for \u}\ > ujc, 



(12) 



which approximates the saturation of scattering rate at 
high frequencies, one obtains A = iaoJc/Tr instead of 0. 

In the next section we show that this problem can be 
partially solved considering a finite frequency range avail- 
able from ARPES experiments. For this range, however, 
the deviation of the bare dispersion from linear becomes 
essential, therefore, we approximate it by a parabola 
e{k) = ujo{l — k'^/kp), which fits well to the band disper- 
sion in given direction obtained within the tight-binding 
model [22|. Here cjo = —/ci?tii?/2 is the bottom of the bare 
band. Within this approximation, the Eqs. I|1U|) and l|ll|) 
can be rewritten as: 
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The fitting machinery is based on Eqs. H13|) . H14|l and 
(121). One can define three steps here. In the two first, 
the real part of the self-energy, for given wq, tOc and n 
(which characterizes the tails, see below), is calculated 
in two ways: (i) Y.'^^^^ by Eq. (jnj); (ii) S'^^ by Eq. ^ 
with subsequent KK transform (|2J). Then, in step (iii), 
the parameters and n are varied until ^^^^.^(cli) 

and I]^^(cj) coincide. In practice, we fit the difference 
'^'disp ^ ^'kk ^'^ ^ small contribution of experimental res- 
olution. 

The calculation of deserves to be considered in 

details. In order to perform a KK transform, high-energy 
tails should be attached to S"(a;) derived from Eq. H14|l . 
Eqs. lO and H12() represent two extremes which can be 
enclosed in a simple analytical expression: 
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as the ultimate cases with — > oo and n — 2 respectively. 
For given n and Wc, we construct Tj"{uj) function in a 
wide frequency range (up to \ujo\ or higher) assuming the 
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particle-hole symmetry: 

E"fw) = ^^w^dthiH) for 1^1 < /-,g^ 

\s:;,„,H for|c.|>c.,„, ^ ^ 

where Um is a " confidence limit" , a maximal experimen- 
tal binding energy to which both the W{u!) and fcm(a;) 
functions can be confidently determined, ^'^idthi'^) '^^l" 
culated from Eqs. (fnji for given wq, ^^^^(w) is fitted to 
'^widthi'^) ii^ the confidence range in order to define a 
and C. Then, we obtain E'^^(w) by KK transform 

Fig- n shows the pairs of and S'(a;) functions 

obtained in such a way for the same loq but for three 
different models; 0, and with n = 4 (dashed, 
dotted, and solid lines respectively). Since KKC = 0, 
in order to simplify numerical calculation, the offset of 
S"(a;) curves is set to S"(wo) = 0- The experimental 
data are taken for an underdoped Bi(Pb)-2212 (UD77, 
Tc = 77 K) at 130 K and the fitting result for which 
tJo = — 0.9 eV is represented by solid-line S'(a;) 0. 

In step (iii), as we mentioned above, the difference 
AE'(a)) = ^Kxi^) - ^'dispi^) should be fitted not to 
zero but to some small but detectable contribution of 
the overall resolution R'{uj) Q. This difference can be 
easily understood by reasoning that finite energy and an- 
gular resolutions mainly effect the MDC's width rather 
than its peak position and that its contribution is fre- 
quency dependent. In order to illustrate this we can take 
into account the overall resolution, _R, as = 
^JYj"{ujY + ■ Then one can consider its frequency de- 
pendent contribution to the imaginary part of as 
the difference between S^-^j^(ijj) and real 

R"{lo) = y/WTrnu^ ^Y."{lo), (17) 

and, due to additivity of the KK-transform, S'^^^^ — 
KKS'^ = KKE" + KKR" = E'^^^^ + R', construct lu- 
dependent contribution to S'^^ as 

R'{uj) = KKR"{uj). (18) 

Although, in principle, the resolution effect R'{uj) can 
be explicitely calculated from known energy and mo- 
mentum resolutions, here we derive it empirically using 
i? as a parameter. It is seen from Fig. ^ that differ- 
ent tails do not affect the energy region < 0.25 eV, 
so, an irreducible difference in the slopes (see Fig. [^J 
A = dYi'jfjf{uj)/duj — dT,'^^gp{Lj) / du) > in the low energy 
range \lu\ < 0.07 eV (while A = at higher energies 0.1 
eV < < 0.2 eV) is a measure of R'{uj). 

In Fig. m we plot R'{uj) setting the offset of to 
zero that gives the value of i? = 0.015 eV. For S"(0) < 
the procedure gives larger R values to accommodate the 
difference in slopes but this does not affect the fact that 
the irreducible difference between S'^^(a;) and 12'^^^j^{uj) 
is caused by the experimental resolution, and depends on 



frequency like it is shown in Fig. 12 it vanishes at zero 
and high frequencies having a maximum around 0.1 eV. 

Thus, we can visualize the fitting procedure as fitting 
the difference AI]'(w) to R'{uj) function. The procedure 
has appeared to be robust with respect to the loq de- 
termination. Fig. 121 illustrates this. First three panels 
show that for a correct value of — —0.9 eV there is 
space for other parameters to vary: different tails can 
be compensated by different ujcS, e.g., for n — 3, A and 
6 in Eq. ^ lOc = 0.34, 0.45 and 0.52 eV respectively. 
On the other hand, at slightly different ojq's (about 10% 
lower and higher, see two right panels), AT,'{oj) cannot 
be fitted to R'{oj) in the whole frequency range. 

OTHER EXAMPLES AND PHYSICAL 
CONSEQUENCES 

In Ref. Q, as mentioned above, we have applied the 
technique to underdoped Bi(Pb)-2212, UD77, in the 
pseudo-gap state (at 130 K). Here we present some other 
examples. 

Fig.Ot shows the results for an overdoped Bi(Pb)-2212, 
OD75, at 130 K. Here, in order to check the correctness 
of the KK numerics, we also plot Y,'^j^(uj) function which 
is obtained by back KK-transform © of I]^jj-(cj). The 
parameters for this fit are: uiq — —0.86 ± 0.03 eV, uJc = 
0.40 ± 0.05 eV, n = 4 ± 0.5. In case of the OD sample, 
the parameters lOc and n are better determined because 
of a higher confidence limit u!,n = 0.45 eV at which one 
can see that S"(ci;) starts to saturate. We note that the 
lower value of uiq compared to the underdoped sample is 
in agreement with the rigid band approximation. 

The demonstrated self-consistency of the analyzed 
data shows, besides the applicability of the self-energy 
approach to superconducting cuprates, that the mea- 
sured spectra belong to a single band and are free of 
influence of any additional features like other bands, su- 
perstructure or A:-dependent background. It has been 
shown recently '3| that although the electronic disper- 
sion along the nodal direction in the bi-layer Bi-2212 is 
not degenerated, i.e. has a finite splitting about 0.05 eV 
for the bare dispersion, the photoemission from the bond- 
ing band is highly suppressed at exactly 27 eV excitation 
energy. At other energies we do not expect that the de- 
scribed fitting procedure will work until the contributions 
of each band can be separated. Fig.|2l) demonstrates this 
showing the " best" fitting result that can be achieved for 
hi' = 38 eV. The difference between AE'((jj) and R'{llj) 
is apparent. 

Fig-Et shows another example of successful application 
of the described procedure to a single-layer Bi(La)-2201, 
OP34, at 40 K. In this case ujq = -0.79 ± 0.05 eV, uj^ = 
0.37 ±0.07 eV, n = 3 ± 0.6. 

The presented examples purpose to illustrate the ap- 
plicability of the self-energy approach to Bi-cuprates. 
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FIG. 3: Real and imaginary parts of the self-energy extracted from the experiment with the described procedure. 



We believe that the described procedure gives a pow- 
erful technique to purify the ARPES data from artifi- 
cial features and to build a strong experimental basis 
for understanding of the nature of electronic interactions 
in cuprates, but still a big work on the data analysis 
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FIG. 4: The resuhs of the fitting procedure for Bi(Pb)-2212 
OD75: (a) real (blue line) and imaginary (red line) parts of 
the self-energy; (b) the experimental (solid red line) and bare 
(dashed line) dispersions on top of the experimentally mea- 
sured quasiparticle spectral weight. 



should be performed. Nevertheless, some conclusions can 
be made even on this stage. 

It is interesting to note that even for UD77 sample, 
for which the saturation of 'S,"{uj) has not been observed, 
it is not possible to reconcile the high-energy behavior of 
S"(cj) with the saturation extreme H12(l [or (jlSfl with n = 
2]. This means that |I]"(cl;)| reaches the maximum and 
starts to decrease at about tOc, and, consequently, 
changes the sign at approximately the same frequency 
(see Fig. For OD and OP samples this conclusion 
is even more strict due to smaller bandwidth. Fig. ^ 
shows the resuhs for Bi(Pb)-2212 OD75: (a) S'(w) and 
I]"(w); (b) fcm(w) and e(fc) on top of the experimentally 
measured quasiparticle spectral weight. 

The fact that ujc is not equal but roughly two times 
less than \ujo\ is consistent with presence of an essential 
electron-electron scattering channel which mainly de- 
termines the lifetime of quasiparticles at high frequencies. 
Other consequences are discussed below. 



MODEL ASSUMPTIONS 

First we focus on two assumptions which have been 
made about the model self-energy: fc-independence and 
particle-hole symmetry. 

It has been mentioned above that the symmetric 
Lorentzian lineshape of the MDC's taken along the nodal 
direction was considered as an experimental evidence 
that the quasi-particle self-energy hardly depends on 
momentum |Q. Recently l^^, it has been noticed 
that the necessary condition for the Lorentzian line- 
shape is dT,"{k,uj)/dk = 0, but dH'/dk can be an w- 
independent constant. This is especially interesting be- 
cause the authors of Ref . ^3 have shown that such a lin- 
ear fc-dependence of S' can explain a non-trivial doping- 
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FIG. 5: Possible particle-hole asymmetry effect on 
(red) and "S^lj) (blue): low-energy (dashed lines, Uc — 0.1 eV) 
and high-energy (solid lines, by Eq. I19|l contributions shown 
on the top of the symmetric self-energy (shaded areas). 



> 

0.05 



0.00 





■ UD77 X ^^^^^ 
O OD75 ^-^-^ 

1 ° °<ft!>°^ 







0.0 -0.1 -0.2 -0.3 

CO (eV) 



FIG. 6: The real parts of the self-energy for UD77 and OD75 
samples at 130 K: solid lines show the result of fitting these 
real parts to 12011 in a frequency range 0.17 eV < < for 
UD77 and 0.12 eV < a; < for OD75; looking down arrows 
mark S'(tj) maxima; the dashed line denotes the 70 meV 
"kink" energy. 



dependent high-energy dispersion observed for a variety 
of cuprates [ig • 

As long as E'(fc,t^) = YJ^{k) + Yl^{ij) and dYJ' jdk = 0, 
/c-dependence of S' does not affect any result of the pre- 
sented analysis except the bare dispersion. In this case, 
the real bare dispersion is just e^'^'^^k) = e{k) — T,'i.{k), or 
^reai _ _ [dT,' / dk) k=k f ■ Although our preliminary 
results, being in agreement with band structure calcu- 
lations ^6] and experimental plasmon dispersion |23j . do 
not support strong fc-dependence of S', it will be interest- 
ing to examine its possible contribution in a wide doping 
range and for different compounds. 

A possible particle-hole asymmetry is another compli- 
cation which can effect the results of the presented anal- 
ysis. In general, one can expect an asymmetry of the 
self-energy due to an asymmetric electron-boson interac- 
tion or as a simple consequence of asymmetric density 
of states. Without considering the origin of the asym- 
metry, we examine its possible influence based on the 
energy scale where it can appear. It is well known that 
because of the possibility to perform ARPES at finite 
temperature one can get the information about quasi- 
particle spectral weight not only below the chemical po- 
tential but also from some region above 01 • For T — 
300 K the MDC width can be measured up to 50 meV 
above Ep, and, within the experimental uncertainty, it 
has appeared to be completely symmetric (e.g., see [U]). 
This means that if there is some asymmetry in the scat- 
tering rate at low energy scale (~ 0.1 eV, a characteristic 
scale which can originate from an electron-boson inter- 
action or from the van Hove singularity in the occupied 
density of states of the hole-doped cuprates), its magni- 
tude is too small to be seen in the < 50 meV energy 
range and, consequently, hardly effects the quasiparticle 
renormalization in the occupied region (w < 0). Fig. [51 
illustrates this: the dashed curves, on top of the sym- 



metric self-energy shown by shaded areas, represent a 
low-energy asymmetric contribution which is too big not 
to be noticed in S"(w) (for < 50 meV) but too small 
to influence 'S'{uj < 0). The solid curves in Fig. |31 present 
the case of high-energy asymmetry that can steam from 
the asymmetry of the bare band We simulate it by 
an asymmetry part in the scattering rate: 

^ \0, uj<0, ' 

where T,'^^^ is determined by ifT^ with coc = 0.45 eV, 
u}c2 — 0.66 eV, n = 4, C = 0. It is seen that although 
the influence of Yi'^(lu) on renormalization at —0.5 eV 
< cij < eV is rather small (can be approximated at 
this stage by a linear contribution with a slope of about 
20% of A) it can be, in principle, detected by more precise 
analysis, in which the influence of the energy and angular 
resolutions is taken into account explicitly. 

PHENOMENOLOGY OF THE KINK 

In Fig. El we plot together the real parts of the self- 
energy for UD77 and OD75 samples at 130 K. Just from 
visual comparison of these data one can conclude that: 
(i) the renormalization for UD77 is considerably higher 
than for OD75; (ii) the energy of the maximum of S'(w) 
for the overdoped sample is lower than for the under- 
doped sample, it is about two times closer to the 70 meV 
" kink" energy; (iii) the kink feature is well defined in the 
underdoped case and becomes weaker with overdoping. 

Following this tendency one can expect that with over- 
doping the 70 meV kink vanishes while the renormaliza- 
tion maximum moves to lower frequencies faking a persis- 
tence of the kink in the whole doping range. Therefore, 
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it is clear that in order to clarify the origin of the kink 
feature a quantitative measure of it is required. 

Keeping the visual definition of the kink as a sharp 
bend of the renormalized dispersion, we formalized it in 
Ref. Q as a peak in the second derivative of S'(a') and 
fitted it to a simple empirical function: 



TT 

X arctan — + arctan 
V 



(20) 



which gives a squared Lorentzian in a second derivative 



K{uj) 



(21) 



Fitting S'(tj) of the underdoped sample in \uj\ < 170 meV 
energy range to this formula we have obtained an energy 
of the kink ujk ~ —63 meV, a kink width [half width at 
quarter maximum of K{uj)] 5 w 30 meV, and a strength of 
the kink / Kduj = AA w 0.65. For the overdoped sample 
« —56 meV, AA ~ 0.45. We believe that a system- 
atic study of this or similar quantitaties as a function of 
doping and temperature will help to find the origin of the 
main electronic interaction in superconducting cuprates. 



CONCLUSIONS 

The self-energy approach is shown to be applicable 
to single band photoemission spectra for underdoped, 
overdoped and optimally doped Bi-based cuprates. The 
demonstrated self-consistency of the procedure opens a 
way to validate the photoemission spectra: the KK-sieve 
can be used to verify the spectra for the absence of the 
band splitting or artificial features. The preliminar anal- 
ysis of the spectra certified in such a way shows that the 
overall renormalization as well as kink in the nodal di- 
rection of Bi-based cuprates is highly doping dependent, 
decreasing with overdoping. In the light of the present 



dilemma about the origin of the main scattering boson 
in the cuprates, a systematic quantitative analysis of the 
KK verified spectra measured at different temperature 
and doping level is indispensable. 
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